Charge transfer in molecular conductors — oxidation or reduction? 
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We discuss the nature of charge transfer in molecular conductors upon connecting to two metallic 
contacts and imposing a voltage bias across them. The sign of the charge transfer (oxidation 
vs. reduction) depends on the position of the metal Fermi energy with respect to the molecular 
levels. In addition, the charge transfer depends on the strength of the coupling (chemisorption 
vs. physisorption) with the contacts. A convenient way to establish the nature and onset of the 
charge transfer and the corresponding features in the FV is to draw an energy level diagram for 
each spin species. Starting from such a level diagram, we argue that transport in the Tour-Reed 
switching molecules, which consist of a central phenyl ring with a nitroamine redox center, involves 
the oxidation of a highest occupied molecular orbital (HOMO)-based level. 



PACS numbers: 05.10.Gg, 05.40.-a, 87.10. +e 

In recent years, there have been many reports 
of conductance characteristics measured for individual 
molecules or small ensembles of oriented molecules [Q. 
Understanding transport in molecular conductors is quite 
challenging because of its special status between solid 
state physics and molecular chemistry. On the one hand, 
the system differs from conventional solid-state devices 
where conduction is understood in terms of alignment 
of a band structure that is largely unaffected by charge 
transfer or the nature of the contacts. On the other hand, 
in contrast to conventional chemistry, the molecule is not 
an isolated unit in equilibrium; coupling with metallic 
contacts makes the system an open one, capable of freely 
exchanging a fractional amount of charge to and from 
the contacts. In addition in the presence of a voltage 
bias across the molecule, the system is out of equilibrium 
and equilibrium statistical mechanical techniques become 
inapplicable. Predicting I-V characteristics of molecular 
conductors involves therefore a suitable combination of 
techniques and insights both from chemistry as well as 
from mesoscopic physics. 

We have recently developed a self-consistent, ab-initio 
procedure for calculating I-V characteristics of molecular 
conductors In this procedure, we self-consistently 
combine the outputs of a standard quantum chemical 
software (GAUSSIAN98) with a non-equilibrium Green's 
function (NEGF) formalism for transport in an open, 
voltage-biased system. The method takes into account 
the hybridization of molecular levels with the contacts 
exactly, and can be used both for weak and strong con- 
tacts. In this paper, we will focus on an intuitive picture 
that can be used to gain qualitative insight into the na- 
ture of conduction through a molecule. This picture is 
particularly useful when the contacts are weakly coupled, 
and issues related to level broadening are minor. The first 
step is to obtain the energy levels of an isolated molecule 



and locate the Fermi energy Ep of the metallic contacts, 
as shown in Fig. ^ (a). Under bias the electrochemi- 
cal potentials /ii^2 of the contacts separate as shown in 
Fig. [^(b) . We wish to address the following questions in 
this paper: (a) which molecular levels are involved in the 
conduction process? (b) Does the molecule get oxidized 
{M M+) or reduced (M M^) at steady state? 

Equilibrium energy level diagram (V = 0). We start 
by obtaining an equilibrium energy level diagram as in 
Fig. 0(a). The energy levels for a given molecule are ob- 
tained using a standard quantum chemical software like 
Hyperchem, MOPAC or GAUSSIAN98. The equihbrium 
Fermi energy of gold is known to be Ep = —5.1 eV. It 
is important to note that semi-empirical programs like 
Extended Hiickel Theory do not give the correct abso- 
lute value for molecular energy levels. In that case, the 
energy levels for gold must also be computed using the 
same program for consistency |^ . 

Which levels conduct?. Fig. |^ shows a schematic 
description of transport through a mesoscopic device. 
At equilibrium (Fig. |l](a)) the metal contacts impose 
their common electrochemical potential on the molecule, 
equal to the metal fermi energy Ep. Coupling of the 
molecule with the contacts broadens each highest oc- 
cupied (HOMO) /lowest unoccupied (LUMO) molecular 
orbital-based level, so that the number of electrons be- 
low the fermi energy is now fractional (the deviation of 
this number from the number of electrons for the iso- 
lated molecule is a measure of charge transfer between 
the molecule and the contacts at equilibrium). When 
a voltage bias is imposed (Fig. |l|(b)), the system is 
driven out of equilibrium and the contact potential splits: 
Ml, 2 = Ep =F eV/2. (Note that in contrast to our earlier 
convention of fixing energy levels and varying the contact 
potentials with bias by different amounts depending on 
the ratio of couplings , we now vary the electrochem- 
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ical potentials as above and let the energy levels adjust 
accordingly). For adiabatic, reflectionless contacts and 
ballistic transport, the contact potentials are imposed 
separately on the right and left-moving electrons As 
long as the two electrochemical potentials lie in the gap, 
all HOMO levels lie below both potentials and are filled 
at zero temperature, the LUMO levels lie above the po- 
tentials and are empty, and there is no current. When 
the bias is sufficient that a level is crossed by one of the 
electrochemical potentials, that level is emptied out (ox- 
idized) from one side but filled (reduced) from the other, 
with a resultant onset of current. Thus only levels which 
lie between the two electrochemical potentials contribute 
to conduction at zero temperature (at finite temperature, 
this range is extended by the thermal energy zbfcsT). 
The onset of current is determined by the proximity of 
the Fermi energy to the molecular level closest to it. 

Oxidized or reduced?. The coupling to the contacts 
of an extended molecule is described by a non-hermitian 
self-energy matrix E calculated in an appropriately cho- 
sen basis This self-energy effectively partitions the 
problem into the molecular part and the contact, and 
leads to a conceptual simplification of the problem. The 
broadenings of the levels (inverse lifetime) are described 
by the anti-hermitian components Ti^2 of Si^2- For fermi 
functions /i_2 corresponding to the two contact potentials 
fii^2, the steady-state occupancy of a molecular level is 
given by / = (/iFi + /2r2)/(ri + r2). If the level hes 
between /ii.2 at zero temperature, /i = and /2 = 1 so 
that the occupancy is given by / = r2/(ri + r2). One 
of the contacts is trying to add electrons (reduce), while 
the other is trying to remove them (oxidize). If the left 
contact is stronger than the right contact (Fi ^ F2), the 
level is emptied (oxidized) while in the opposite case the 
level stays filled (reduced) . In either current flows 

due to the competition between oxidation and reduction. 
Typically for a self-assembled monolayer, one end of the 
molecule (usually consisting of a thiol group chemisorbed 
on a gold surface) has a strong contact, while the other 
end is physisorbed and has a weak contact. This means 
that for a positive bias on the physisorbed end, the 
molecule remains neutral, while for opposite bias, charg- 
ing occurs from the strongly contacted chemisorbed end 

For a high work function metal, Ep usually lies near 
the HOMO level. This is an issue that clearly needs a lot 
more attention, since the precise location of Ep depends 
on the model for the contacts, as well as the method 
of calculating the molecular Hamiltonian. A third ter- 
minal (gate) [|j can help resolve some of these issues. 
One of the consequences of a self-consistent evaluation 
of energy levels is that as we raise the fermi energy, the 
levels float up by an amount equal to the charging en- 
ergy (for a solid, this is roughly given by the capacitance 
of a sphere of radius equal to the size of the electronic 
wave function). For a molecule with strong coupling to 



at least one contact, or for a molecule in solution, the 
effective electronic size is infinite, so that the charging 
energy is small, and classical STM theory holds. For a 
weakly coupled system, charging effects lead to the lev- 
els following the Fermi energy, so that the conduction 
tends to remain HOMO-based. We verified this for the 
phenyl-dithiol molecule (PDT) where we couple the out- 
puts of a GAUSSIAN98 [|| evaluation of energy levels 
using a Becke-3 parameter exchange and Perdew-Wang 
91 correlation (B3PW91) approximation [ p^pT| , with a 
non-equilibrium Green's function (NEGF) evaluation of 
the density matrix |^,^ , and iterate to convergence self- 
consistently The basis set used was the Los Alamos 
National Lab set for effective core potentials of the double 
C type (LANL2DZ) ||l|,|ll. The converged energy lev- 
els increase with increasing Ep at the rate of 4eV/Volt, 
which corresponds to the capacitative charging energy of 
a sphere of size 5A. 

Level diagram under bias (V ^ 0). As an illustra- 
tive example, we have plotted the energy levels of the 
Tour-Reed molecule 2'-amino- 4,4'-di[ethynyl]phenyl-5'- 
nitro-l-[thioacetyl] benzene containing the nitroamine re- 
dox center iQ, as a function of voltage applied to the 
sulphur end (Fig. ^). We consider an extended molecule 
including three surface gold atoms bonded to sulphur and 
obtain energy levels analogous to Fig. |[ including the 
effect of voltage-induced Stark shifts of the molecular 
orbitals. Exchange and correlation effects are included 
by performing a B3PW9I calculation with a LANL2DZ 
basis set using the GAUSSIAN 98 software. In the 
Tour-Reed experiment, one end of the molecule is self- 
assembled onto gold using sulphur "alligator clips" , lead- 
ing to a strongly coupled chemisorbed contact, while the 
other end is a carbon atom weakly coupled to an evap- 
orated gold contact. Charging occurs when the electro- 
chemical potential /ii of the sulphur end (dotted line in 
Fig. H) crosses a molecular level. From Fig. ||, we see 
that for positive bias on sulphur, we first cross a HOMO 
level by the stronger sulphur contact, which oxidizes the 
molecule (the first LUMO level is actually a nonconduct- 
ing gold-based level, as we discuss below). For negative 
bias, the HOMO level is crossed by the weaker carbon 
contact, so the level stays filled and there is no charging. 
This means that the I-V should be strongly asymmetric, 
as is indeed observed. Q The sulphur contact potential 
Hi crosses the HOMO level at around 2 V, whereupon 



*For a spatially symmetric molecule, such contact-induced 
charging asymmetries alone determine the asymmetry in the 
I-V llq]. However for spatially asymmetric molecules such 
as the Tour-Reed example here, the I-V asymmetry is much 
larger, and given predominantly by the unequal Stark shift- 
ing of the intrinsic molecular wavefunctions (the asymmetric 
Stark shifting is seen in Fig. Q). 



2 



the molecule is oxidized from M to M+. This is im- 
portant, since the precise nature of the charging process 
is an important question raised in the literature Jl6| , p^ . 
Although cyclic voltammetric curves seem to suggest re- 
duction processes as being operative, one has to remem- 
ber that in such measurements counterions are present 
in solution, and substantially compromise any charging 
effects that may dominate in solid state. 

When M changes to M+, we need a whole new set of 
energy levels (Fig. |^). Two changes happen: (i) since 
one electron of a particular spin is pulled out by the con- 
tact, the spin degeneracy is removed; (ii) secondly, the 
removal of one electron of given spin leads to charging 
effects which cause all the levels of obtained from 
GAUSSIAN98 to float down by about 2 volts. This is 
unphysical, because GAUSSIAN98 models an isolated 
molecule with image charges at infinite distance. In the 
presence of image charges on the contacts, the charging 
energy will be lower. We therefore raise the energy levels 
of M"*" so as to align the energy level of the semi-occupied 
molecular orbital (SOMO) with the HOMO of the neu- 
tral species {M). Such an alignment is motivated by 
our assumption that in the absence of correlations and 
molecular reorganization, the ionization potential of M 
equals the electron affinity of Af+. This assumption is a 
point that requires further work, and should come auto- 
matically out of a proper ab-initio theory that includes 
the effect of image charges on the contacts. Note that 
in order to obtain spin-dependent charging effects that 
shift all levels except the SOMO substantially, one needs 
to do a spin-unrestricted (UHF/LSDA) calculation that 
breaks the spin degeneracy of the energy levels. 

It is interesting to note that for an extended molecule, 
there are some gold-based levels. These can be identi- 
fied either from their wave functions that are localized 
strongly on gold, or from their energy levels; the levels 
run roughly parallel to the electrochemical potential /ii 
of the contact on the chemisorbed side (dotted line in 
Fig. H). For small molecules these gold levels, as well 
as metal- induced gap states (MIGS) that arise when a 
self-energy is included in our calculation, end up con- 
tributing to the conductance through direct tunneling 
processes. For large molecules such states hybridize with 
molecular wave functions, but don't contribute directly 
to tunneling, so they can be ignored. However the gold 
levels can still influence the positions of the other levels 
through level-repulsion, so the contacts affect the voltage 
at which a HOMO level is crossed and oxidation occurs. 

We have shown that in self-assembled monolayers such 
as the Tour- Reed molecules containing the nitroamine re- 
dox center, charge transfer occurs due to oxidation from 
the strongly contacted end involving the level closest to 
the Fermi energy. We are currently investigating the 
mechanism of switching thereafter. Understanding the 
behavior of a molecular switch is important from the 
point of view of molecular electronics. Several switching 



molecules have been reported [ liljiq l, and various pos- 
sible mechanisms suggested ||r7,|l9|. A charged species 
could deform in an applied field, and the electromechani- 
cal effect could switch off the current. If the added charge 
is localized on a sidechain of the molecule, it could rotate 
the central benzene ring attached to it in the presence of 
the field, distortiing the conjugated electron wave func- 
tion enough to switch current off. A more complicated co- 
operative effect involving the whole self-assembled mono- 
layer could also arguably lead to a current switch. Prior 
to a full self-consistent ab-initio calculation of transport 
properties including contact effects, it is important to 
develop an intuitive understanding of transport in the 
isolated species. An energy level diagram as illustrated 
above can serve as a quick intuitive tool to answer some 
basic questions such as the onset of current, whether an 
electron is being added or removed, or the existence of 
an asymmetry in the I-V. 
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FIG. 1. Simplified description of conduction in a meso- 
scopic device: (a) Left: At equilibrium (V = 0) the metal 
contact electrochemical potential is given by /ii,2 = Ef- The 
molecular levels are discrete, but end up being broadened by 
coupling with contacts, (b) Right: Under bias (V> 0) the 
system is out of equilibrium as the contact potentials split, 
/ii,2 = Ef T e.V/2. Only the levels lying between ^1 and ^2 
are filled from one side (oxidized) and emptied from the other 
(reduced) with a resultant conduction of current. Levels lying 
below or above both electrochemical potentials stay filled or 
empty and do not contribute to current at zero temperature. 
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FIG. 2. Top: Tour-Reed molecule with nitroamine 
side-group, extended to include three gold atoms from the 
contact surface bonded directly to sulphur. The molecule 
has a strong chemisorbed bond at the sulphur side and a 
weak physisorbed one at the other. Bottom: Energy levels 
for the neutral molecule (Af) as a function of voltage ap- 
plied to the sulphur end. The levels Stark shift with applied 
voltage. There are some gold-based levels running roughly 
parallel to the electrochemical potential on sulphur /ii (dot- 
ted line), that do not contribute to conduction. The Fermi 
energy Ef ~ —5.1 eV is closer to the HOMO levels than the 
LUMO ones (ignoring the trivial gold-based levels). At pos- 
itive bias, the HOMO level is crossed first with the strong 
sulphur contact (dotted line) around 2 V, which empties the 
level out (oxidation). Subsequently, we draw the energy lev- 
els for the oxidized species M"*", which are spin nondegenerate 
(red and black lines). At negative bias, the HOMO levels are 
crossed with the physisorbed carbon contact (dashed line), 
which does not charge up the system due to the weakness 
of the contact coupling. The resulting I-V is expected to 
be strongly asymmetric, and exhibits charging effects arising 
from an oxidation step at around 2 V. 
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